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ABSTRACT 

We present a first analysis of the g-mode oscillation spectrum for the white dwarf (WD) primary of GW Lib, 
a faint cataclysmic variable (CV). Stable periodicities have been observed from this WD for a number of years, 
but their interpretation as stellar pulsations has been hampered by a lack of theoretical models appropriate to 
an accreting WD. Using the results of Townsley and Bildsten, we construct accreting models for the observed 
effective temperature and approximate mass of the WD in GW Lib. We compute g-mode frequencies for a range 
of accreted layer masses, M acc , and long term accretion rates, (M). If we assume that the observed oscillations 
are from I = 1 g-modes, then the observed periods are matched when M ss 1 .02M Q , M acc « 0.3 1 x 10~ 4 Mq and 
(M) « 7.3 x 10 -1 1 Mq yr" 1 . Much more sensitive observations are needed to discover more modes, after which 
we will be able to more accurately measure these parameters and constrain or measure the WD's rotation rate. 

Subject headings: binaries: close — novae, cataclysmic variables- stars: dwarf novae — white dwarfs 



1. INTRODUCTION 

Dwarf Novae (DN) are the subset of CVs with low time- 
averaged accretion rates (M) < 10~ 9 Mq yr" 1 and thermally 
unstable accretion disks that lead to sudden accretion events 
which interrupt the otherwise quiescent state. In quiescence, 
the UV (and sometimes optical) emission from the binary is 
dominated by light from the WD surface, allowing for a mea- 
surement of the WD's r e ff. These T e gs are much hotter than 
expected for a WD of the age of the bin ary (a few Gyr) and 
must be related to accretion iSionlll99^) . Calculations of the 
heatin g of the deep interior of the WD by the prolonged accre- 
tion ( Townslev & Bildsten 2004) explains the observed values 
of r e ff and yields a unique relationship betwe en (M), the WD 
mass, M, and the orbital period P or b ( Townslev & Bildsten 
2003). 

GW Lib is o ne of the shortest o rbital period CVs known, 
Porb = 77 min (Thorstensen et al. 2002 ), and therefore has 
a very low (M) ~5x 1O -11 M vr" 1 ^Townslev & B ildsten 
2003). Only one disk outburst has been observed from GW 
Lib JGonzalea ll 9 8 3) . Much later photometric observations 
dur ing quiescence led to the discovery of periodic variabil- 
ity ( Ivan Zvl et alJ2000l) similar to that of isolated W Ds which 
pulsate due to non-radial g-modes (see lBradlevlll998l for a 
review of the DAV WDs). The highest signal to noise pho- 
tometric observations comprise two weeks of single-site data 
taken in 1998 with some supplement from other longitudes 



TABLE 1 
Observed Periods 



System 


Principal Period, P 




(seconds) 


GWLib 


236, 377, 646 


SDSS 1610 


345, 607 


SDSS0131 


~ 330, ~ 600 


SDSS 2205 


~ 330, ~ 600 
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( Ivan Zvl et al.120 04). There are three clear periodicities, listed 
in Tabled with some evidence for mild period variability. 

Three addition al CV WD pulsators have b een found 
JWarner & Woudd 120031 IWoudt & Warned 12004) by taking 
photometric time series of DN identified in the Sloan Digi- 
tal Sky Survey (SDSS; Szko dv et al.ll2003l) which show WD 
spectral features in the optical. These objects and their dom- 
inant periods are also listed in Table ^ At this discovery 
fraction we expect ~15 CV WD pulsators total by the end 
of the SDSS. Seismology of the CV WD pulsators provides 
the prospect of well-determined masses for these systems and 
the parameters we derive for GW Lib are the first step in this 
process. 

Mass determinations are especially interesting due to pos- 
sible implications for progenitor systems of Type la super- 
novae. Though t to be accr eting WDs near the Chandrasekhar 
mass (iHillebrandt & Niem ever 2000), these progenitors are 
closely related to the CV population, but the precise nature 
of this relationship is unknown. Important clues in this mys- 
tery lie in how the masses of CV primary WDs change over 
the accreting lifetime of the binary, but progress is hampered 
by the difficulty of measuring C V primary masses ( P atterson! 
1998). This mass evolution, as well as probing the M acc di- 
rectly with seismology, is also important for determining how 
much, if any, of the original WD material has been ejected 
into the ISM in classical nova out bursts, contributing to the 
ISM metallicity JGehrz et alll998l) . 

We begin in §2 by discussing the observed properties of the 
WD in GW Lib and constructing accreting WD models that 
are consistent with the observations. Section 3 discusses the 
g-mode properties for the accreting model of GW Lib, and 
in Section 4 we show what can be learned about the accreting 
WD from the observed mode periods. We conclude in Section 
5 with a discussion of future work. 

2. OBSERVED PROPERTIES OF THE WHITE DWARF IN GW LIB 

Though the periods make it clear that the oscillations in GW 
Lib are non-radial g-modes, without secure identification of 
the radial and angular quantum numbers (n and £), it is diffi- 
cult to carry out the seismology. Hence, we start by using the 
binary's observed properties to constrain a few of the WD's 
properties, thus reducing the range of possibilities we need to 
consider in our modeling efforts. 
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FIG. 1. — Parameter space constraints from non-seismological observa- 
tions. Distance, UV flux and UV spectral measurements constrain M (diag- 
onal shading) and gravitational wave radiation provides a (mass dependent) 
lower limit on (M) (solid line). Actual models considered here are indicated 
by the circles. The * indicates our best fit to the observed pulsations. 



2.1. Limits from Distance and T e ff 

We use three observations: (1) the parallax (2) a r e ff-logg 
relationship determined from the UV spectrum, and (3) the 
P 01 b. The parallax of GW Lib has been measured from the 
ground ( Thorst ensenl 120031) and gives a distance of 104^0 
pc. Current UV spectroscopic observations can only con- 
s trai n and logg to a well-determined linear relationship 
(see Howell et al. 2002 for a discussion). In the case of GW 
Lib, r eff /K = 14700 + 2000 log w g& where g$ is the surface 
gravity measured in 10 8 cm s" 2 (Szkodv et al. 2002). By us- 
ing this, the UV flux measured from the same observation, the 
distance, and a WD mass-radius relation we find M = 1 .03 to 
1.36M© from the range of allowed distances. This constraint 
excludes the diagonally shaded regions in Figure 

Another parameter needed is the time-averaged accretion 
rate, (M). For each M, P OI b sets the orbital separation, and 
with a mass-radius relat ion for the Roche Lobe filling com- 
panion (we use Kolb & Baraffe's 1999 for a low mass main 
sequence star) this allows a derivation of a mimimum (M) 
due only to gravitational radiation from the orbit. This limit 
is shown by the solid line in Figure [2 and can be high if the 
system has passed through the period minimum, so that the 
companion is a sub-stellar object which is out of thermal equi- 
librium and has a lower me an density. 

Optical spectroscop y (ISzkodv. Desai. & Hoard! 120001 
iThorstensen et al.| 120 02) has given lower r e ff (11000 K and 
13220 K respectively), however both of these studies used 
gravities much lower than the allowed range derived above 
(logg = 8 and 7.4 as opposed to 8.6 for M = 1.0M Q ), and 
were unable to effectively account for contamination by the 
quiescent accretion disk, both of which would lead to low 
fitted r e ff. The UV spectrum is the most reliable indicator of 

?eff. 



2.2. The Accreting White Dwarf Structure 

The WD interior structure is from iTownslev* & Bildsten 
(2004, hereafter TB). We chose a simple compositional struc- 
ture consisting of a solar composition accreted layer of mass 
M acc on the WD core, an equal mixture of 12 C and 16 0. For 
calculation of the buoyancy properties, a smooth transition 
region of 0.2 times the local pressure scale height was put in 
between these layers (see TB for a discussion of the diffusion 



timescales). With this compositional structure, a WD model 
is parameterized by M, T core , (M) and M acc , and compression 
of material by accretion powers a surface luminosity L. TB 
related r core to (M) by finding the equilibrium state where be- 
tween classical novae outbursts (as M acc is growing) the WD 
core would suffer no net heating or cooling. For a fixed T e ff at 
the observed value, TB also relate (M) andM acc , with a higher 
M acc implying a lower (M) for the same T e ff. Using these two 
constraints we are left with M and M acc as free parameters, 
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where we only consider M a , 
nova ignition mass. 

This smaller grid of models are shown by the circles in 
Figure Q] and have M/M© = 1.0,1.05,1.1,1.15, and 1.20, 
each with M acc = 0.1, 0.3, 0.5, 0.7, and 0.9M ign . The igni- 
tion masses for these models were M; gn = 1.40, 1.32, 1.22, 
1.14, and 1.08 x 10 _4 M Q . The model at M = 1.05M o and 
M acc = 0.3M ign = 0.40 x lO^M© has T eS = 16070 K (that im- 
plied for this M by the UV spectrum), (M) = 6.4 x 10" U M Q 
yr" 1 , and r coie = 5.9 x 10 6 K. The downward trend in (M) with 
increasing M is due to L decreasing with the WD radius with 
r e ff constrained to the measured relation. 

3 . G-MODES IN GW LIB ' S WHITE DWARF 

In a nonrotating star, all variables can be decomposed 
into spherical harmonics Yi, n (9,(f)). We ignore the perturbed 
gravitational potential (the Cowling approximation) and work 
in the adiabatic approximation. The linearized momentum 
and energy equations are then written in terms of the ra- 
dial displacement £,-(r) and the Eulerian pressure perturbation 
Sp(r) = pip(r) as (See e.g. eq.14.2 and 14.3 in lUnno et alJ 
1989) 

N 2 , , 

:—^-(N 2 -LU 2 )C r , 
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Here u> = 2tt/P is the mode frequency, P is the pe- 
riod, c = (T\p/ p) x l 2 is the adiabatic sound speed, g = 
GM(r)/r 2 is the downward gravitational acceleration, and 
Af 2 = — g (d\n p/ dr d\n p / dr) is the square of the Brunt- 
Vaisala frequency. For boundary conditions, we impose zero 
Lagrangian pressure perturbation, Ap/p = ip-g^ r = 0, at the 
top of the model and £ r = at the interface between the solid 
core and surrounding liquid. This interface occurs at the 
freezing point for a Coulomb solid (Bi ldsten & CuflerJl 19951 
Montgome ry & Wingeti n"999) and g-modes cannot penetrate 
into the solid core as their frequency is too low to excite sig- 
nificant shearing motion of the Coulomb solid. A solid core 
is present in the GW Lib model due to the high mass. 

The propagati on cavity f or the observed short wavelength 
g-modes ( Unno et al. 1989) in GW Lib is bounded from above 
by lo 2 = L 2 = 1(1+ l)c 2 /r 2 , the Lamb frequency, and by the 
solid core below. The propagation diagram for the model 
with M = 1.05M© and M acc = 0.3M ig n = 0.40 x lO^M© is 
shown in the upper panel of Figure |5] The large peak in 
N at log 10 p ~ 19 is due to the change in mean molecular 
weight in the transition layer from the solar composition ac- 
creted envelope to the C/O core. In the roughly constant 
flux envelope, A^ 2 ~ g/z where z is the depth from the sur- 
face. In the degenerate core, A^ 2 ~ (g/H p )(kbT /Ep), where 
Ef is the electron Fermi energy and H p = p/pg is the pres- 
sure scale height. In the propagation zone of the wave, the 
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FIG. 2. — Comparison of propagation diagram (top panel) and WKB inte- 
grand (bottom panel) for an accreting model (solid lines) and non-accreting 
model (dashed line) with the same composition and T c ff. The vertical solid 
(dashed) line on the right indicates the liquid-solid boundary. The solid hori- 
zontal lines are the observed frequencies and the dotted horizontal line is the 
breakup frequency. The accreting model has M = 1.O5M0, M acc = 0.3Mi„ n = 
0.40 X 1O^M , T eff = 16070 K (that implied for this M by the UV spec- 
trum), (AT) =6.4 X 1O~ U M yr _1 ,and T corc = 5.9 X 10 6 K. The non-accretor 
has r corc = 9.3 X 10 6 K. 



WKB dispersion relation for low frequency g-modes gives 
a radial wavenumber k r = (N/ui)[l(l+ l)] 1//2 /r. The quan- 
tized WKB phase ip accrued by the wave in a given region is 
(p = J drk r =[l{l+\)] x l 2 uj- 1 jNdr/r so that uj oc jNH^dlnp. 
The integrand is shown in the bottom panel of Figurey] rep- 
resenting the number of nodes per decade in pressure. Shown 
are two WD models: an accreting model with M = 1.05M Q 
and M acc = 0.3M ig „ = 0.40 x 10 _4 M Q used in our mode analy- 
sis (solid line), and a non-accreting model with the same T e g 
(dashed line) and composition. 

Not shown in Figure [2] is the difference due to the enve- 
lope mean molecular weight, /i, from the case of a pure H 
or He envelope to one of solar compositon. Such a change 
is reflected in the periods as P„ oc /i 5 , since the envelope is 
well-approximated by an n = 4 polytrope. The non-accreting 
model also has a higher T core than the accreting model by 
about 50%, leading to two important effects: (1) the WKB 
integrand has a higher value in the core for the non-accreting 

1/2 

model (uj oc r core ), and (2) the solid core is smaller, pushing 
the inner boundary condition deeper into the star. Both effects 
directly influence the observed periods, and period spacings, 
hence it is essential to use a WD model including the effects 
of compressional heating and nuclear burning, rather than a 
passively cooling WD model. 



4. INFERRING THE WD PROPERTIES FROM THE PERIODS 

The existing o ptical time-series photometry of GW Lib 
llvan Zvl et al.l l2004) consists of 7 time series taken in 1997, 
1998, and 2001. Here we focus on the best of these, the two 
weeks of data from May 1998. Table [2 represents our esti- 
m ate of the three perio ds. From the O — C phase plots shown 
bv lvan Zvl et al.l (1200 4). the periodicity near 646 s varies with 
v = oj/2tt ~ -10 -11 Hz s -1 , and the phase is tracked by the ob- 
servations. For that near 377 s, the phase is lost, and therefore 
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FIG. 3. — Contours in the root mean squared difference between the 
mode periods observed and the closest modes in the model, AP lms . The 
best solution lies in the lower left hand corner where M = 1 .02Mq , M acc = 
0.23Mjg„ = 0.31 X 10~ 4 Mq. The dot-dashed line demonstrates approximately 
how M acc /Mj gn scales with M for a fixed period and n. 



it is inconclusive whether there are multiple closely spaced 
periods or just unresolved variability. The highest frequency 
is quite stable, and though close to a sum is unlikely to be a 
combination frequency. Splittings ~ 1 /iHz are ubiquitous 
in the Fourier transforms, but their origin is unclear. The 
Doppler shift due to the orbit should have a magnitude of 0.5 
to 1 .3 /iHz for the 646 and 236 s periods respectively. The fre- 
quency resolution set by the inverse length of the time series 

is also ~ /iHz. 

As mentioned by van Zvl et alJ (|2004), the mode frequen- 
cies might drift due to the cooling of the material accreted in 
the last DN outburst. The magnitude and sign of this drift can 
be approximated using the cooling time of the freshly accreted 
material and its fractional contribution to the WKB integral 
discussed above. For a recurrence time of 20 years (the time 
since the last outburst) with (M) ~5x 10" n M Q yr" 1 , AM ~ 
10" 9 M Q was deposited in the last outburst. For M = l.OM 
the base of this layer has ph — 2x 10 14 dyne cm" 2 , and a ther- 
mal time of r th = c p T b AM/L ~ 1 yr (p fe /10 14 dyne cm" 2 ) 12 , 
implying that this layer relaxes to close to the static solution 
in a few years after the outburst, though a small amount of 
cooling continues. As can be seen from the lower panel in 
Figure|2l the mode frequencies are largely determined deeper 
in the WD; however, this outer layer does have a contribu- 



tion J Ndz oc pi ° oc T^, for an n = 4 polytrope. This ac- 
counts for about 20% of the whole integral, and leads to 
v ~ 0.2vf/2T = -v/lOm ~ -10" 12 Hz s" 1 for v = 1/646 s, 
close enough to the observed drift to warrant further calcula- 
tions of this effect, which we do not undertake here. Shorter 
period (lower n) modes reside deeper in the star, making them 
less vulnerable to transient heating and cooling effects. 

We have calculated the g-mode spectrum of model WDs 
with parameters on the grid indicated in Figure^for I = 1 . By 
interpolating within this grid we compute the model periods 
P„(M,M acc ). The root mean square period difference between 
the observed three modes and the models is defined as AP 2 ms = 

^^ =1 rrun{„).[(P n -f > i ) 2 ]/3, where the observed periods have 
been indexed with z = 1,2,3. The contours of AP rms are shown 
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in Figure [3] implying a best fit solution with AP lms = 1.8s 
forM= 1.02M© and M acc = 0.23M ign = 0.31 x 1O" 4 M . The 
mode identifications for this best fit model are n = 3, 8, and 
17 for the three observed periods. These do not appear to 
correspond to any particular mode trapping pattern reflected 
in the mode kinetic energy. For this model, we can predict 
additional mode periods yet to be observed; up to n= 17 these 
are 141, 191, 234, 268, 289, 310, 351, 377, 400, 432, 463, 
492, 519, 553, 586, 617, 647, 678 s. 

The contours for Af rms do not close around a single so- 
lution, implying additional data is needed to get better con- 
straints. There are two other areas in the M-M acc plane which 
also provide fairly good matches (P lms < 3 s) to the observed 
periods. These are centered near (M/M©, M acc /Mj gn ) of 
(1.06, 0.53) and (1.15, 0.6), and extend in the M acc direction. 
Each of the three minima represent a different set of mode 
IDs, the shallower ones having n = 4, 9, 19 and 5, 1 1, 22 for 
the 236, 377, and 646 s modes. The M, M acc degeneracy can 
be simply explained. For a mode trapped in the envelope, 
P„ ocn/ J eav Ndz oc nM~^! c l M~ 23 for a poly trope of index 4 and 
an assumed scaling R oc M" 1,8 , suitable forM = 1.0 to 1.2M . 
Thus fixing n and P n defines the relation M acc /Mi gn oc M" 215 , 
where M; gn oc M" 1 5 , leading to the degeneracy shown by the 
dot-dashed line in Figure|3] 

5. CONCLUSIONS AND UNCERTAINTIES 

We have found non-rotating accreting WD models which 
fit the periodic variability observed in GW Lib, implying 
M=1.O2M , = 0.31 x 1O _4 M and (M) =7.3 x 10- n M Q 
yr" 1 . This model is quite simple, but since it is capable of ac- 
counting for the observed periods, a more complex model is 
not called for at this time. A number of important issues were 
not dealt with, such as the presence of a residual He layer due 
to inefficient nova mass ejection, or our treatment of the ac- 



creted layer-C/O boundary. More detailed non-adiabatic cal- 
culations will be undertaken to address mode excitation. The 
most essential omission in our calculations is the WD rotation, 
since it is expected to be spun up by accretion. If the WD is 
rapidly rotating, the model parameters derived here are likely 
to be inaccurate. Better observations are necessary for this 
next essential parameter to be constrained. These observa- 
tions should accomplish two goals (1) to clearly characterize 
the substructure or variability of the 377 s and 646 s periods, 
(2) probe for lower amplitude signals. A larger number of 
identified modes are essential for conclusive seismology. 

Rotation may significantly modify the mode frequencies. 
Since the maximum rotation rate for our GW Lib mod- 
els is ~ 1 rad s" 1 and the observed mode frequencies are 
ui ~ 10~ 2 rads _1 , spin frequencies above 1% breakup will 
cause large frequen cy shifts for observed modes. In fact 
ISzkodv et alJ (12002) estimate vsini < 300 kms" 1 , allowing 
spin rates up to 6% of breakup even for high i. Rotation 
rates this large lead to a qualitative change in mode prop- 
erties. When the spin frequency O sp i n < 0.5o>, rotation can 
be treated as a small perturbation. Also, observationally this 
rotational frequency splitting allows a determination of the I 
quantum number. If, however, i! sp j n > 0.5w, the Coriolis force 
has a "nonperturbative" effect on the mode frequencies, mod- 
ifying the relati on between the frequency a nd quantum num- 
bers (Bildsten. Ushomirskv. & Cutler 1996). Hence accreting 
WDs may provide a unique proving ground for techniques of 
seismology for rapidly rotating stars. 
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